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Abstract. One of the most spectacular phenomena in physics in terms of dynamical 
range is the glass transition and the associated slowing down of flow and relaxation 
with decreasing temperature. That it occurs in many different liquids seems to call for 
a "universal" theory. In this article, we review one such theoretical approach which 
is based on the concept of " frustration" . Frustration in this context describes an 
incompatibility between extension of the locally preferred order in a liquid and tiling 
of the whole space. We provide a critical assessment of what has been achieved within 
this approach and we discuss the relation with other theories of the glass transition. 



1. Introduction 

One of the disappointing aspects of the study of the glass transition is that after many 
years of scholarly effort there remains wide divergence in the basic ways in which relevant 
phenomena are envisaged. By "basic ways to envisage phenomena" we mean the physical 
models upon which relevant theories of structure and dynamics can be built, i.e., the 
zeroth-order descriptions which with but minor elaboration should enable one to describe 
universal or species-independent phenomena over a wide range of control variables such 
as temperature and pressure. Perhaps such great differences of viewpoints exist because 
no such zeroth-order description can be constructed, but we do not believe this to be the 
case. What might be described as the " anomalous slowing down" of relaxation and flow 
in liquids with the approach to the glass transition appears to be a sufficiently general 
phenomenon to be describable by a "universal" theory. 

We review in this article one such theoretical approach which has been developed 
over the past two decades: it is based on the idea that the physical mechanism which 
is responsible for glass formation is the ubiquitous presence of "frustration" in liquids. 
"Frustration" in this context describes an incompatibility between extension of the local 
order preferred in a liquid and tiling of the whole space. 
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After briefly summarizing the phenomenology of glassforming liquids in light of 
what we think to be the relevant physics (section |2J, we introduce the concept of 
frustration. We first describe what is probably the most studied example, that of 
frustrated icosahedral order in metallic glassformers, and we summarize the curved- 
space approach of geometric frustration (section EJ). In section |U we discuss the 
statistical mechanical description of frustration in liquids, with both the underlying 
assumptions and the possible candidates for providing a minimal theoretical model at 
a coarse-grained, mesoscopic level. We next present the evidence and the consequences 
of a generic property shared by the frustrated models for the glass transition, that of 
"avoided criticality" (section EJ). In section® we review a phenomenological scaling 
approach to supercooled liquids built about an avoided critical point, the frustration- 
limited domain theory, and in section we summarize the results of extensive computer 
simulation studies of 3-dimensional Coulomb frustrated lattice models. The following 
section is devoted to discussing additional phenomena brought about by frustration, 
the emergence of a complex free-energy landscape characterized by a proliferation of 
long-lived metastable states and the presence of topological excitations (defects) whose 
effective dynamics may be strongly constrained; this allows us to make contact with 
other theoretical approaches of glass formation (section |HJ). Finally, section |U] presents 
some concluding remarks. 

This article is intended to provide a critical assessment of what has been achieved 
within the frustration-based approach of glassforming liquids, discussing the limitations 
and the open questions. 

2. Phenomenology: what is there to be explained? 

What is meant by the "glass transition" is not always well specified, but usually it 
refers to the passage from liquid-like behavior to behavior characteristic of amorphous 
solids. The "glass transition" so defined is not a thermodynamic, nor even a dynamic, 
transition; it is a point of dynamic arrest arbitrarily but narrowly specified by the 
available experimental techniques for studying flow or relaxation. We focus here on the 
approach to this glass transition, when the liquid or the polymer can still be considered at 
equilibrium (although it may be " supercooled" and therefore in metastable equilibrium, 
the most stable phase being a crystal). 

The distinctive property of glassforming liquids and polymers is the dramatic 
slowing down of relaxation and flow with decreasing temperature P] El HI] ■ The dynamical 
range of this phenomenon is probably one of the most spectacular in physics: viscosity 
and main (a) relaxation time may increase by some 15 orders of magnitude for a 
mere 30% decrease in temperature. Especially striking is the fact that for most liquids 
(notable exceptions are network forming systems with strong directional bonding such 
as silica glasses), the slowing down occurs much more rapidly that one would have 
anticipated by extrapolating the behavior of the "ordinary" liquid, say above its melting 
temperature T m . For instance, extrapolation of the high-temperature data for the 
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Figure 1. Super Arrhenius temperature dependence of the viscosity of liquid ortho- 
terphenyl: logarithm (base 10) of rj versus inverse temperature 1/T. The straight line 
is the extrapolation of the behavior in the "ordinary" liquid range (above T m ); whereas 
the actual T g is around 243-fT, the extrapolated one is at a much lower temperature 
around 50 — 70K. (Data taken from references cited in 0]). 



viscosity of liquid ortho-terphenyl, which are well described by a simple Arrhenius 
dependence, r\ ps r] 0C exp[E 00 /(kBT)], would imply a glass transition temperature 
(conventionally defined as the point at which the viscosity takes a given value of 10 13 
Poise) at very low temperature, some 150 to 200K below the actual T g » 243K: 
see figure [0 This points towards the existence of a crossover temperature T* below 
which the dynamical properties of a liquid change from an "ordinary", Arrhenius- 
like temperature dependence to an "anomalous", stronger than Arrhenius (or super- 
Arrhenius) one. Such an anomalous dependence, which is observed in all kinds of 
liquids and polymers irrespective of their detailed molecular properties [U El Ej, suggests 
a collective phenomenon involving cooperative motion of a large number of molecules. 

The notion of "fragility", introduced by Angell to classify glassforming liquids 
and polymers [H], characterizes the degree of super- Arrhenius behavior, the more fragile 
the glassformer the greater the super-Arrhenius character. The strong temperature 
dependence of the viscosity and a-relaxation time of glassforming liquids is often 
described by a Vogel-Fulcher-Tammann formula, n « ^exp [DT /(T — To)], D being 
a measure of fragility. This description implies the presence of a dynamical singularity 
at a temperature T lower than T g . (However, T is always found quite far below 
T g : at least 40 K for the fragile liquid ortho-terphenyl, see figure [TJ) An alternative 
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Figure 2. Crossover from Arrhenius to super Arrhcnius behavior for several 
glassforming liquids: effective activation (free) energy E(T) versus 1/T. Both E(T) 
and T are divided by the crossover temperature T* . (Data taken from references cited 
in Refs.^ and 

way of representing the data is to focus on the effective activation (free) energy, 
E(T) = k bT In (ri/rjoo). In this representation, the crossover behavior shows up quite 
clearly, as illustrated in figure El for several liquids. 

Another characteristic of the approach to the glass transition, which becomes more 
manifest in the temperature range below the crossover discussed above, is the spatially 
"heterogeneous" nature of the relaxation. The most easily detected, yet not uniquely 
interpret able, signature of such heterogeneity is a stretching of the relaxation functions, 
i. e., a deviation from simple exponential dependence at long times (or the equivalent 
broadening in frequency space for the dynamic susceptibilities); this is often represented 
by a stretched exponential in time, (f)(t) ~ exp(— (t/r Q ) /3 ) where r a is the a-relaxation 
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time. More direct experimental evidence for the presence of supermolecular dynamically 
correlated regions in supercooled liquids has been provided in the past ten years (see the 
reviews [HE! El)- This, too, suggests collective behavior; however, the measured (linear) 
size of the correlated regions never exceeds 5 to 10 molecular diameters at T s [10j. 

An additional puzzling feature of the glass transition problem is that the 
spectacular change in the dynamical properties does not come with a related growth in 
thermodynamic correlations 121 E]. No convincing evidence has been so far reported 
for a significant increase in structural correlations: the evolution of the static pair 
structure factor appears weak and featureless even in very fragile glassformers. The only 
noticeable observation is the decrease of the excess entropy, defined as the difference 
between the entropy of a supercooled liquid and that of the associated crystal, with 
decreasing temperature ^l] . This decrease appears correlated with the slowing down of 
the relaxation, being larger for more fragile glassformers [ID] 

We can then summarize our view of the phenomenology of supercooled liquids as 
follows: glass formation seems to be a collective or cooperative (we use the two terms 
interchangeably) phenomenon, which would account for the rather universal trends 
observed in many different systems, but cooperativity occurs on a limited, mesoscopic 
spatial scale. 

As we have mentioned in the Introduction, the absence of a unique reading of the 
salient aspects of the phenomenology of glassforming liquids goes with the absence 
of a widely accepted theory of the glass transition. There is indeed a variety of 
theoretical approaches which often appear at odds with one another. The one feature 
that most theories have in common, though, is the idea that the slowing down of 
relaxation and flow in supercooled liquids (especially in the fragile ones) is spectacular 
enough and general enough to be describable by a universal theory and that, as is 
quite natural in physics, universality is associated with the presence of one or several 
underlying critical points: the putative critical points may be dynamic [121 EH1 EH EE] 
or thermodynamic [ini E3 EH EH EH El 122]; but, as no divergences or singularities 
are actually observed, those critical points are postulated to be either unreachable 
[Tol E3 EH EH EH EE1 I2H 122] i- e., occurring below the experimental T g , or avoided 
[121 ED], being only present in an ideal system "near" to the real liquid (see figure EJ). 

In this article, we review one such theoretical approach, which is based on the 
concept of frustration. Frustration, meant here to describe an incompatibility between 
the locally preferred order in a liquid and the global requirements for tiling of space, 
is an appealing candidate for generating collective (cooperative) behavior on a limited 
spatial scale and bring in static spatial correlations only at a multi-particle level that 
cannot be probed by usual experimental techniques. 
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Figure 3. Schematic representation of (critical) theories of the glass transition in 
a virtual multidimensional diagram: the physical axes (temperature, volume) are 
complemented by abstract dimensions (frustration, hopping mechanisms). Avoided 
critical points: T* = T c , critical point in zero frustration; T^ ICT , dynamic singularity 
of mode coupling theory. Unreacheable critical points below T g : Tk, entropy crisis 
or random first-order critical point; T = 0, dynamic critical point for kinctically 
constrained models; point J, jamming transition. (See text for references.) 

3. Geometric frustration in simple atomic systems: the curved-space 
approach 

The concept of frustration was introduced by Toulouse in 1977(211 in the context of spin 
models to describe situations in which one cannot minimize the energy of the system 
by merely minimizing all local interactions. A simple illustration is provided by an 
Ising spin model on a triangular lattice with antiferromagnetic interactions between 
nearest-neighbor spins: an elementary "plaquette" of the lattice made by one triangle 
of nearest-neighbor sites is frustrated because the three antiferromagnetic interactions 
along the bonds forming the plaquette can never be satisfied simultaneously (see figure 
13} . Spin glasses are a well-known example of frustrated systems; however, in this case 
frustration is induced by the presence of quenched disorder, due for instance to frozen-in 
impurities[23]- This inhomogeneous and externally imposed frustration is not relevant 
for supercooled liquids in which glassiness and heterogeneous behavior are self-generated. 
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Figure 4. Elementary plaquette of a triangular lattice. When interactions between 
nearest neighbors Ising spins are antiferromagnetic, one can never satisfy (i.e. anti- 
align spins) more than 2 bonds at the same time. 

To emphasize the difference, frustration in the context of liquids has been given several 
qualifying adjectives, such as "uniform", "geometric", "topological", or "structural". In 
what follows however we shall simply refer in most cases to "frustration" without adding 
any epithet. 

What has often been taken as the paradigm of (geometric) frustration in amorphous 
systems is the case of icosahedral order in simple one-component liquids in which the 
atoms interact through spherically symmetric pair potentials. The ground state of four 
atoms is a perfect tetrahedron, with the atoms sitting at the vertices, and twenty such 
tetrahedra can be combined to form a regular 13-atom icosahedron (figure EJ). Frank [2*K] 
was the first to stress the importance of local icosahedral order in simple atomic liquids 
and liquid metals. He showed that the most stable cluster made of a central atom and 
a shell of 12 neighbors is indeed an icosahedron, and not an arrangement associated 
with the actual crystalline phases, bec, fee, or hep. In 3-dimensional Euclidean space, 
however, this locally preferred structure (tetrahedral or icosahedral) cannot propagate 
freely to tile the whole space: this is what is meant by geometric frustration. The global 
ground state of the system is instead an fee or hep crystal. 

The concept of frustration is more easily grasped by constrasting the case of 
spherical particles in 3 dimensions with other situations [26, 27]. Systems of spherical 
particles in 2-dimensional Euclidean space, i. e. disks on a plane, are not subject to 
frustration: the locally preferred structure is a regular hexagon, with one atom at the 
center and 6 neighbors at the vertices, and this structure can be periodically repeated to 
form a triangular lattice. Similarly, aligned hard cubes in 3-dimensional Euclidean space 
form simple cubic arrangements that are both locally and globally preferred, therefore 



The frustration-based approach of supercooled liquids and the glass transition 



8 




Figure 5. Manifestations of frustration for tetrahedral/icosahedral order in 3-d atomic 
systems: (a) the 5-fold rotational symmetry of the icosahedron shown here is not 
compatible with translational periodicity, (b) The distance d between the atoms of 
the first shell is slightly larger than the distance a between the central atom and its 
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excluding any frustration effect 

As illustrated in figure frustration can be envisaged from various different angles: 
as a global constraint, e. g., the incompatibility between the 5-fold rotational symmetry 
of the icosahedron and translational periodicity (figure OJi), or, more in line with the 
usual definition of frustration (see above), as the impossibility to let all the atoms 
simultaneously sit in positions corresponding to the minima of the pair potentials 
between nearest neighbors (figure 03). 

An important observation is that frustration can be turned off by "curving" 
space j2H EH- One can play with the metrics and topology of the underlying space in 
order to let the local order propagate freely and generate a regular tiling or "tessellation" . 
This can be understood for instance by considering figures Et> and c: if space can be 
curved so that the defect angle made by 5 tetrahedra packed along a bond as well 
as the mismatch between the two interatomic distances d and a in an icosahedron 
vanish, propagation of tetrahedral/icosahedral order becomes possible. Indeed, a perfect 
icosahedral phase can be formed on the surface S 3 of a 4-dimensional hypersphere with a 
radius tuned to be 5/rr = 1.5915... times the preferred interatomic geodesic distance (in 
curved space). This ideal icosahedral structure is called the polytope {3, 3, 5}[27| I30j. 

The core of the curved-space approach of amorphous (liquid and glassy) phases, 
developed on one hand by Sadoc, Kleman, and Mosserij23 EH EH EH ESI an d on the 
other by Nelson and coworkers [HH ESI ESI EE1 EH EH] and SethnajSHl HOj, consists in 
using this ideal structure as a reference state (for reviews see Refs. (2SJ EH EH 132] • Going 
back from the reference state to the actual configurations in "flat" (Euclidean) space 
necessarily forces in topological defects that perturb the ideal order. All sorts of defects 
can be generated, but it was argued that the most relevant for the physics of atomic 
glasses are "disclinations" [321 EH EH ES] : disclinations are associated with the breaking 
of rotational, here bond-orientational, symmetry (their counterpart for translational 
symmetry are the "dislocations"), and they are line defects in 3 dimensions [13]. Such 
topological defects are expected to form a disordered network, the actual amorphous 
phases being composed of locally ideally ordered regions interrupted by defects [44J. 
Interestingly, the network of defects (disclinations) may itself become ordered at low 
temperature, thereby leading to "defect-ordered phases" whose best examples in the 
case of metallic systems are the Frank-Kasper phases (23 E3 HE]- 

Whereas Sadoc and Mosseri focused on the most efficient way to produce dense 
amorphous packings by decurving the polytope {3, 3, 5}, Nelson and coworkers |3(Il EH 
EE1EH] developed a statistical mechanical approach of frustrated icosahedral order. They 
built a set of local order parameters Q^(x) by projecting a local particle configuration 
centered on position x onto the surface of a tangent 4-dimensional hypersphere with 
appropriate radius to accommodate the polytope {3, 3, 5} and, much in the spirit of the 
Ginzburg-Landau theory of freezing, they derived a free-energy functional of this local 
order parameter [35]. 

F[Q t ] = [ d?x{\z, | Z^Q,(x) | 2 +\ n | Q,(x) | 2 +0(Q/)}, (1) 
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where the index \i is associated with the spatial coordinates in ordinary 3-D 
Euclidean space (e. g., // = x, y, z for cartesian coordinates); the local order parameter 
Q^(x) transforms according to the representation of dimension (£ + l) 2 of SO (4), Zg 
and re are phenomenological parameters, and 0(Q/) denotes cubic and higher-order 
invariants built with Q^. The most relevant value of I for describing icosahedral order 
is 12, but t = 20, 24, 30, 32, ... are also allowed by symmetry (£ — corresponds to the 
mean particle density). 

Frustration is introduced via the so-called "covariant derivative" entering in the 
gradient term. Following a suggestion made by Sethna [231101, it is defined by requiring 
that the gradient term be minimum when two nearby local particle configurations are 
related by "rolling without slipping" the tangent hypersphere carrying the reference 
icosahedral template along the path joining the configurations. This leads to[36 

Z^Q,(x) = ^Q,(x) - «cL$Q*(x), (2) 

where the Lq2's denote the generators of the £0(4) rotations in the (0, /i) plane (0 
denotes the direction in the 4-dimensional space embedding £ 3 which is perpendicular 
to the tangent Euclidean space described by the cartesian coordinates fj, — x,y, z) and k 
is the inverse radius of the hypersphere £ 3 compatible with the polytope {3, 3, 5}. Here 
and in the following, the Einstein convention is used, in which one sums over indices 
whenever they are repeated. The covariant derivative is frustrated because it cannot be 
minimized (i. e., set equal to zero) everywhere. 

Despite its many appealing features, this curved-space approach of frustrated 
icosahedral order has led in practice to meager results concerning the glass transition 
phenomenon itself. Several reasons may explain this shortcoming. First, the models, 
such as the £0(4) uniformly frustrated theory described above, are rather complex 
and have never been fully analyzed: for instance, the topological defects have been 
characterized |3S1 OH dEj and the static structure factor of metallic glasses has been 
reproduced [H3 EE], but no attempts have been made to study the dynamics. 

Secondly, difficulties may arise from the fact that local icosahedral order in simple 
atomic systems appears strongly frustrated. The typical distance between defects, 
expected to be of the order or less than the radius of the reference hypersphere K~ l , 
is small: one or two particle diameters at most, which precludes a separation between 
a mesoscopic, coarse-grained description of collective effects and a detailed account of 
the microscopic properties. (This may also explain the controversy about the extent of 
icosahedral order and associated correlations in amorphous phases composed of spherical 
particles.) It is worth noting that simple atomic systems are not good glassformers and 
appear quite "non-fragile" with an essentially Arrhenius temperature dependence of 
the viscosity and a-relaxation timejUj; this is also true for the Lennard- Jones models 
studied in computer simulations which, in spite of the weakness of the interatomic 
bonding, are much less fragile than typical molecular glassforming liquids 48J. 

Finally, if glass formation is indeed a "universal" phenomenon whose most salient 
properties are largely independent of molecular details (in a sense discussed in the 
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preceding section), icosahedral order is likely to be too material-specific. Possible 
generalizations to describe frustration in liquids are reviewed in the next section. 

4. Statistical mechanics of frustration in liquids 

Application of the concept of frustration to liquids boils down to three propositions 
which are plausible and, as discussed above, reasonably well established in the case of 
systems of spherically symmetric particles, but still remain at the level of postulates 
in the case of molecular liquids, mixtures, and polymers. Those propositions can be 
summarized as follows: 

1) A liquid is characterized by a locally preferred structure (LPS) which is different 
than that of the crystalline phases. This LPS is an arrangement of molecules that 
minimizes some local free energy (an operational definition of such a local free energy 
is given in Ref. |49j). much in the same way as the local tetrahedral/icosahedral order 
in liquids of spherical particles. Description of such a LPS requires the knowledge 
of multi-particle densities, beyond the usual one-body density, and of multi-particle 
correlations, beyond the usual pair correlations. For instance, even in simple atomic 
liquids, information on local bond-orientational order involves at least 2-body densities 
and 4-body correlations [5(3 EH E2|- For this reason, the signature of a LPS is hard to 
detect in common diffraction experiments in which one has only access to the (pair) 
structure factor of a liquid or a glass. 

The LPS may change at high enough pressure: this is the case for instance of 
water and tetrahedrally bonded liquids for which the low-pressure LPS that tends to 
minimize the local energy by satisfying the maximum number of directional (hydrogen- 
or covalent) bonds goes over to a different structure at high pressure where packing 
effects dominate over bonding effects [SHI El E3- However, at a given pressure or 
density, the LPS does not change with temperature; it just becomes more and more 
favored as temperature decreases. 

2) The LPS characteristic of a given liquid cannot tile the whole space. This 
incompatibility between local order and global space filling is precisely the content of 
the concept of frustration. Were it not for frustration, the LPS could propagate through 
the whole space and form an "ideal" ordered phase that would supersede the actual 
crystalline phases. Instead of this, the actual crystal is formed via a strong first-order 
transition that requires a rearrangement of the local structures, a cost that is more than 
compensated by the free-energy gain due to the global tiling of space. So, as suggested 
by Frank [2*5], the presence of frustration is the possible physical ingredient allowing 
supercooling of a liquid at temperatures below the melting point. 

3) It is possible to construct an abstract reference system in which the effect of 
frustration is turned off. At a microscopic level, frustration implies that the atoms or 
the molecules cannot simultaneously sit in the minima of all pairwise interactions with 
their nearest neighbors. By allowing one to modify the topology and the metrics of the 
underlying space, one can get rid of frustration. This is the essence of the curved-space 
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approach of icosahedral order (see above). In an abstract parameter space, one can now 
think of frustration as a tunable control variable. However, a given liquid is of course 
characterized by a given level of frustration, a given LPS, and a given reference space. 

Note also that in the absence of frustration, "ideal" ordering is likely to proceed via 
a continuous or weakly first-order transition, quite differently than the standard melting 
in 3 — D flat space. For instance, the melting transition of parallel cubes, a system with 
no frustration between local order and global tiling in ordinary 3-D space, has been 
shown to be continuous [28 . Similarly, Nelson and his coworkers [3T)j I3*T| EHJ EE] have 
argued that bond-orientational ordering is likely to occur via a continuous or weakly 
first-order transition in 3 dimensions, even in the presence of cubic invariants in the 
Landau description of the svstem jo^ loT]. (Icosahedral ordering in the polytope {3, 3, 5} 
on S 3 [E3 does not occur by a true phase transition at all since the underlying space is 
the surface of a 4-D sphere of finite radius.) 

The next step is to build a statistical mechanical description of glassforming liquids 
based on the above propositions. We will take here the point of view that if indeed, as 
argued in discussing the phenomenology, the viscous slowing down of supercooled liquids 
and the resulting glass transition are universal, progress can be made by relying on a 
coarse-grained, mesoscopic approach that incorporates the effect of frustration while, at 
least in a first step, overlooking a detailed description of the LPS, of the reference space, 
and of the precise mechanism by which frustration operates |58|. 

A minimal theoretical model of frustration in glassforming liquids should satisfy, 
on top of the overall consistency with the physical picture developed above, some 
prerequisite properties, such as (i) frustration should be a tunable parameter, which 
allows one to investigate the potential connection between frustration and "fragility" 
in glassformers, and (ii) frustration should be uniform, i.e., the same at every point 
of space and, consequently, not due to the presence of quenched disorder. We review 
below some of the paths that have been, and still are, followed. (Note that in all these 
descriptions, the actual crystal is excluded from the picture, in order to focus on the 
liquid and supercooled liquid phases.) 

4-1. Coupling to a non-Abelian gauge background 

The statistical mechanical treatment of icosahedral order put forward by Nelson and 
his coworkers (see above) can be extended by ignoring the precise reference of the local 
order parameter (and of the LPS) to the ideal tiling of the polytope {3, 3, 5} and that 
of the frustration strength to the radius of curvature of the hypersphere S 3 compatible 
with this ideal tiling. 

In a continuum, field-theoretical description, the generic form of the free-energy 
functional reads \59 1 

F[Q] = J d?x{ l -Z | D„Q(x) | 2 + l -r | Q(x) | 2 +0(Q 3 )}, (3) 

where the index fi is associated with the spatial coordinates in ordinary 3 — 



The frustration-based approach of supercooled liquids and the glass transition 



13 



D Euclidean space; the local order parameter Q(x) transforms according to the 
representation of dimension n of a given non-Abelian continuous Lie group, Z and r 
are phenomenological parameters, and 0(Q 3 ) denotes cubic and higher-order invariants 
built with Q. 

The "covariant derivative" entering in the gradient term is defined as 



considered non-Abelian group. The A^s are coefficients depending on both "spatial" 
indices fi and "internal" indices a; they are independent of position and play the role 
of a frozen, uniform gauge field (see below); the coupling strength g characterizes 
the magnitude of the frustration. (We recall that the Einstein convention is used for 
summing indices.) 

The thermodynamic and structural properties of the system can been obtained 
from the partition function 



where (3 = l/iksT), whereas studying the relaxation properties requires the 
introduction of some local dynamical rules that reflect, at a coarse-grained level, the 
underlying motion of the atoms or molecules. 

For example, in the theory developed by Nelson and summarized in section El 
the non-Abelian group is SO (4), the main local order parameter transforms according 
to the representation of 5*0(4) associated with angular momentum / = 12, i. e., 
of dimension n = (I + l) 2 = 169, a denotes the 6 rotation planes in 4-D space, 
(0,x),(0,y),(0,z),(x,y),(x,z),(y,z); A° is equal to one when a corresponds to the 
planes (0,/i) with \i = x,y,z and zero otherwise, and g is equal to k*, the inverse 
of the radius of the "ideal" reference hypersphere. A related, but somewhat simpler 
approach due to SethnajS3 HO] corresponds to the I = 1 (i. e., n = 4) representation 
of SO (4) (see also Ref. [SI])- Lattice versions of the above free-energy functional have 
also been proposed (BOIES], which could make possible computer simulation studies of 
such theories. 

In the present description, frustration is embodied in the specific form of the 
covariant derivative. When the coupling strength g is set to zero, the reference, 
unfrustrated model is chosen such that there is a continuous or weakly first-order 
transition to a symmetry-broken phase in which the order parameter takes a non-zero 
value (this phase describes the ideal order formed by extending the LPS to the whole 
space). When g is different from zero, frustration arises from the fact that the covariant 
derivative cannot be made zero everywhere. For instance, if one builds configurations 
by forcing the covariant derivative to vanish along straight lines starting in all directions 
from a given position x, one finds that the covariant derivative differs from zero along 
any closed circuit encircling x; this leads to a large free-energy cost associated with a 



Z^Q(x) = ^Q(x) - i(?A«L a Q(x) 



(4) 



where the L a 's denote the generators of the representation of dimension n of the 




(5) 
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non-zero gradient term. (Actually, as discussed in more detail below, this cost grows 
super-extensively with the linear size of the region around the arbitrary chosen point x, 
up to some intrinsic frustration length set up by the value of g.) 

In the framework of differential geometry |61| I2T| H3J making the covariant 
derivative vanish along a curve amounts to parallel transporting the local order 
parameter (which is a vector or a tensor) along this curve. In the curved-space approach 
of Sethna and Nelson, this is done by "rolling" a tangent 4-D sphere S 3 with the ideal 
icosahedral template along the chosen curve in Euclidean space. Frustration then means 
that one cannot extend such a parallel transport to the whole space. The reason is that 
the covariant derivative has a non-zero curvature; this curvature is directly related to 
the magnitude of the difference between the initial and final values of the local order 
parameter after parallel transport along a closed circuit. The presence of a non-zero 
curvature (not to be confused with the curvature of the physical Euclidean space which 
is of course zero: in the Sethna-Nelson picture, the non-zero curvature of the covariant 
derivative is precisely equal to that of the reference 4-D sphere) comes from the non- 
commutativity of the generators L a , i. e., from the non-Abelian property of the gauge 
background. 

An important consequence of frustration is the necessary presence of defects: the 
ideal order that forms the ground state of the unfrustrated system (g = 0) can no 
longer propagate freely through the whole space, and low-temperature configurations 
must contain an irreducible density of defects (e. g., disclination and dislocation lines, 
grain boundaries, etc.). Such "topological" defects can be studied by means of the 
homotopy theory [o^j. as done by Nelson and Widomj^n] for their 5*0(4) theory of local 
icosahedral order. 

For completeness, it is worth interpreting the above class of models in the language 
of gauge field theory [4*^1 142j. The local order parameter Q(x) represents the matter field 
that is "minimally coupled" to the non-Abelian gauge field. In the present case, the 
gauge field is frozen and uniform (hence the term "gauge background" used to describe 
such a theory). The associated physical field, defined as the so-called covariant curl 
F^y = — i [D^ D v ] , is then given by F^ u = ig 2 A^A^[L a , Lp\. It is non-zero because 
of the non-commutativity property of the L a 's associated with the non-Abelian nature 
of the "gauge group" under consideration. Contact can be made with the formalism 
of differential geometry and tensor calculus by noticing that identifies with the 
curvature tensor |4*3]. Note that the present theory is only globally invariant under 
transformations of this group. If necessary, a full blown gauge field theory in which 
local invariance is now enforced via a fluctuating, space-dependent gauge field could be 
developed; attempts in this direction have been made EH EHJ EH| • 

4-2. Uniformly frustrated spin models 

In developing their approach based on geometric frustration for simple liquids and 
glasses, Nelson and his coworkers [33 1213 EH] proposed an analogy with the effect of 
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an applied, uniform magnetic field in extreme type-II superconductors [U7|. In such 
systems, their is always at least partial penetration of the magnetic field, no matter how 
small the field: this latter induces an irreducible density of defects (flux lines or vortices) 
of the same "sign" . The ground state is then an Abrikosov lattice in which those vortices 
form a periodic array, in a similar way as Frank-Kasper phases in which disclination lines 
form a periodic network are the expected ground states of simple metallic systems. The 
curvature mismatch that generates frustration in the curved-space approach is replaced 
here by the applied magnetic field. 

The analogy can be made somewhat more precise in two dimensions. As discussed 
in section El there is no frustration in that case for liquids formed by one-component 
particles interacting via spherically symmetric pair potentials: the hexagonal LPS 
can propagate through the whole space to form a triangular lattice and the freezing 
transition is either weakly first-order or splits into two continuous transitions separated 
by an intermediate hexatic phase characterized by bond-orientational order but no 
translational order. Geometric frustration can now be introduced by curving space, i e., 
by placing the liquid on a hyperbolic plane H 2 , a surface of constant negative curvature 
— /t 2 [26 , EH1 ESj • Let focus on bond-orientational order. Frustration then comes from the 
fact that if one measures the bond angle formed by the vector joining two neighboring 
atoms with respect to a reference axis, this reference axis changes when it is parallel 
transported on the curved surface. (On a curved manifold, there is no global notion of 
parallelism, so that comparison of reference frames at different points in space requires a 
rule for parallel transport of vectors along curves.) A non-zero curvature implies that the 
reference axis does not go back to its original value when transported along any closed 
circuit. In the long wavelength, continuum approach (also called "hydrodynamic" in this 
context), the Ginzburg-Landau free-energy functional describing local hexatic order on 
a hyperbolic plane can be written as a modified gradient term[35j, 

F[n] = l -Z H J d 2 xyfgjxj | D,n"(x) | 2 , (6) 

where n(x) is a unit vector tangent to a "bond" centered at point x, g(x) is the 
determinant of the metric tensor (^(x) that is appropriate for the coordinate system 
chosen for the hyperbolic plane, Zh is a phenomenological hexatic stiffness constant, 
and the covariant derivative is given by 

D/(x) = 5/(x) + r;(x)n A (x), (7) 

where the 's are the so-called connection components that define parallel 
transport on H 2 ; here they are simply the Christoffel symbols obtained from derivatives 
of the metric tensor and associated with the Levi-Civita connection. In a local cartesian 
coordinate system centered at an (arbitrary) origin, one has (^(x) = 5^ — K 2 x fl x u + 
0(k 4 x 4 ) and T^(x) = T^ A (x) = —k 2 5\^x u + 0(k 4 x 3 ) for distances much less than the 
frustration, or curvature, scale, re -1 . 
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Consider now the Ginzburg-Landau free-energy functional for an extreme type-II 
superconducting film in a perpendicular magnetic field B[67\: 

/h 2 lei 1 

d 2 x{^ I (0„ - —A^m*) | 2 +-r | if)fa) | 2 +0(^ 4 )}, (8) 

where if) fa) is the complex superconducting order parameter and A(x) is the vector 
potential whose curl, d A A = B, is a fixed constant (where B = Bz is perpendicular to 
the plane (x, y)). 

Rewriting if) fa) in terms of two real fields, if) fa) = if) 1 fa) + iif> 2 fa), and dropping 
terms other than the (modified) gradient one, Eq.(JHJ) can be recast as[35j 

m = £ / ^ 1 ( + r ^( x )^ A ( x )) I 2 ' ( 9 ) 

with the "connection" now defined as fa) = (2ei/tvc)A l]L fa)e\i l , e\ u being 
the usual antisymmetric tensor with e 12 = — e 2 i = 1 and v4^(x) being equal to 
A^fa) = — (l/2)Be fia x rT , due to the uniform property of the magnetic field B (see 
above). Comparing Eqs. ()6l7|) with Eq.(JH} and using the expansions of g^ufa) and T^fa) 
for distances less than the frustration scale, one sees a great similarity between the two 
systems, with B playing the role of k 2 ; the mapping however is not exact because the 
tensorial content of the various objects is not quite the same (in particular the gauge 
group U(l) associated with superconductors is Abelian, which is not the case for that 
associated with the connection on H 2 ) and because curvature introduces an intrinsic 
scale k~ 1 beyond which the correspondence breaks down. 

The analogy between the two kinds of systems is also striking when looking at a 
dual picture in terms of defects. The relevant defects in the case of hexatic ordering 
are disclinations, which are points in two dimensions. When frustration is weak, i. e., 
for small curvature, and at low temperature, it has been shown that local hexatic order 
on a curved surface (here, the hyperbolic plane H 2 ) can be described in the continuum, 
"hydrodynamic" limit by a free-energy functional for the defect (disclination) density 
s(x),PHZ01 

F[s) = NE core + J d 2 x^g~fa) j d 2 y^g~fa)(sfa) + K 2 )Gfa, y)(s(y) + k 2 ), (10) 

when the total number of defects is equal to N, E core is the energy penalty 
associated with creating one defect, and the density of defects is defined as sfa) = 
fa/S) Yu%=i n Qi^fa — x i)' w ith Qi — -F2, ••■ the topological charge of the defects; 
Gfa, y) is the inverse Laplacian on the hyperbolic plane, which goes as the coulombic 
interaction ln(r/ao), with r the distance between two defects and a$ the core size of the 
defects, up to a screening length played by the intrinsic frustration or curvature scale 
The above expression is obtained by requiring "charge neutrality", which means 
that j d 2 x^f gfa) (sfa) + k 2 ) = or else that the mean defect density is equal to the 
Gaussian curvature — k 2 up to some trivial prefactor. This charge neutrality expresses 
the frustration effect of curvature that forces an irreducible density of defects of the same 
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sign (here, point disclinations with a negative topological charge). As stated previously, 
this is quite similar to the frustration effect of an applied uniform magnetic field on 
a superconducting film, in which an irreducible density of vortices of the same sign is 
induced. Actually, at low temperature and in the continuum limit, the dual description 
of Eq.fjHJ in terms of defects leads to the following free-energy functional [7T ] I72~ ] 173] : 

F[s] = NE core + l -Z J d 2 x J d 2 y(s(x) - /)G(x,y)( S (y) - /), (11) 

where G(x, y) is the inverse Laplacian in flat space, i. e., goes as the coulombic 
interaction log(\ x — y | /a ) at large distances, and / oc B 2 . Here too, charge neutrality 
imposes that the mean density of vortices is equal to /. 

We have made the above rather long digression to illustrate Nelson's analogy 
between glassforming liquids and uniformly frustrated systems. Taking this analogy for 
granted, one may envisage as possible minimal theoretical models for glass formation 
simple uniformly frustrated modelsJTI], either in the continuum version described above 
or in a lattice version, which may be more convenient for computer simulation studies. 
The simplest of such models is the two-dimensional uniformly frustrated XY model, 
whose Hamiltonian is 

H = -jJ2 cos ( e i - °j ~ A a) ( 12 ) 

(ij) 

where the sum is over distinct pairs of nearest-neighbor sites, 0j is the angle of the XY 
spin at site i, and Ay is the line integral of the vector potential along the bond (ij); the 
sum of the A^-'s taken in a clockwise direction over the bonds forming an elementary 
plaquette of the lattice is fixed to be J2n^ij = ^ n f f° r eacn plaquette (/ G [0, 1/2]). 
The frustration parameter / is simply related to the square of the amplitude of the 
uniform applied field (see above). By duality transformation [75] . one may alternatively 
focus on the associated defect (vortex) description: this corresponds to a 2d Coulomb 
lattice gas in a uniform background charge — /, whose Hamiltonian is just the lattice 
version of Eq. (jllj) . 

Simulations of the 2d uniformly frustrated XY model and of the associated Coulomb 
lattice gas model have been performed for a large value of frustration, /* = (3 — v / 5)/2 
corresponding to the so-called "full frustration", and glassy behavior reminiscent of 
supercooled liquids has indeed been observed [THJ [77]. 

In three dimensions, one can generalize the Coulomb lattice gas to a lattice model of 
interacting vortex loops, in which vortex lines are forced by frustration in all directions 
of space (as is required to model a glassforming liquid, but no longer corresponds to an 
extreme type-II superconductor in an applied magnetic field[78j): 

H=^ Yl (^-/Mlxi-Xjl)^-/) (13) 

where Sj M is the integer vorticity on bond /i emanating from site i, which forms 
continuous lines (loops or infinite lines), and G(\xi — Xj|) is the lattice 3d Coulomb 
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Figure 6. "Squashing" of an ideally ordered region of size L from curved space to flat 
space; I denotes the intrinsic curvature (i. e., frustration) scale. 

interaction that goes as l/|xj — Xj| at large separations. Charge neutrality imposes that 

At this point one should stress that whereas structural and thermodynamic 
properties are simply obtained by studying the partition function computed by 
integrating over the defect degrees of freedom, the dynamics is trickier: the mapping 
from the local order parameter, which evolves with simple local dynamics, to the 
associated defects may imply quite complex kinetic rules for the defects (see section IHH . 

4-3. Competition between effective interactions: (screened) Coulomb 
frustrated models 

Frustration describes a competition between a local tendency to order and opposing 
constraints. It is then tempting to describe this situation at a coarse-grained level by 
means of competing effective interactions acting on different length scales. The situation 
seems quite clear for the ordering tendency: the corresponding interaction is short- 
ranged and the associated local order parameter characterizes the liquid LPS. Although 
things are not as clear as far as the frustrating effective interaction is concerned, there is 
a good rationale for taking this interaction as long-ranged. This is what we argue now. 

A robust feature of the effect of frustration is the fact that, at low temperature, 
forcing the extension of the LPS (i. e., the ideal order) over a region of size L induces 
a super-extensive free-energy cost, a strain energy in a continuum elastic picture, that 
grows as L d+2 in dimension d. A heuristic and somewhat simple-minded, argument can 
be put forward as follows. Imagine a reference curved space with intrinsic curvature 
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scale £ which is (ideally) ordered at low temperature (see figure |HJ). What is the effect of 
"squashing" an ordered region of size L, with L much larger than the typical molecular 
size but smaller than £ (this implies that frustration is small), onto flat, Euclidean space? 
the free-energy density of this region in flat space can be expanded as 



The first term of the above expression corresponds to the free-energy gain associated 
with ideal ordering in the absence of frustration ($oW < 0); the second term is zero, 
$i(£) = 0, because one can choose the flat region locally tangent to the curved space. 
The last term describes the dominant contribution to the cost generated by frustration. 
As a result, the corresponding contribution to the free energy of the region of volume L d 
in flat space goes as ($2(£)/£ 2 )L 2+d for a space of dimension d. This "super-extensive" 
cost prevents the extension of the LPS in Euclidean space. 

The argument can be checked on the models described in the two above subsections. 
The 2-dimensional case is straightforward. Consider the free-energy functional described 
in Eqs. fl6l7|) . In this case, ideal order consists in a hexatic phase which is present in the 
unfrustrated Euclidean space and corresponds to a uniform order parameter n(x) = n . 
The free energy of a region of size L <C £ = k~ x with n(x) = n on the hyperbolic plane 
is then given at leading order in k£ by 



where we have used the fact that at leading order, « 5^ and (x) = 
— k 2 5\(j,x u . Here too, and up to possible logarithmic corrections, there is a free-energy 
cost that grows as L d+2 with d = 2. The same reasoning applies to the extreme type-II 
superconducting film in a perpendicular magnetic field. 

In terms of the associated defect picture, with defects being either point 
disclinations or point vortices, the above result can be simply interpreted: forcing the 
ideal (hexatic or superconducting) order in a region of size L when frustration is present 
amounts to forcing the defect density to be zero in that region; due to the coulombic 
nature of the interactions in Eqs. (fTUjl and (fTTj) , this then leads to a super-extensive cost 
that goes as k 4 L 4 or f 2 L 4 (plus logarithmic corrections). A main difference between 
hexatic order on the hyperbolic plane and superconducting order in a uniform magnetic 
field is the presence of a screening effect in the former case at distances larger than the 
intrinsic frustration (curvature) scale. 

The 3-dimensional case of models with a coupling to a non-Abelian gauge 
background is slightly more involved. Consider for illustrative purpose the 5*0(4) theory 
of frustrated icosahedral order with the local order parameter Q^(x) with £ = 12 (in the 
following we drop the index £). Ideal order now exists in curved space. Following the 
interpretation of Sethna 40^ and Nelson- Widom . particle configurations in Euclidean 
space that most closely mimic the ideal order are constructed by "rolling" a reference 



$(L, = $ (Z) + $!(0(y) + $ 2 (/)(y) 2 + 0( A 3 ). 



(14) 




(15) 
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polytope {3, 3, 5} between adjacent points. It is then convenient to introduce a change 
of variable 



Q(x) = e^^Q(x) (16) 

so that the order is now measured relative to a reference polytope that has been 
rolled in straight lines emanating from the origin in all directions. One can then try to 
force the relation Q(x) = Q in a region of linear size L around the origin, with Q 
minimizing the potential part of the free-energy functional (i. e., \t | Q(x) | 2 +0(Q 3 )). 
The free-energy cost is associated with the gradient term J d 3 x | D At Q(x) | 2 . 

By appying the Baker-Campbell-Hausdorff formula 

e A e B = e A + B +hi A ^+T2(i A 'iAB]]+lB,[B,A}])+- ( 17 ) 



to 



where we have made explicit the sums over repeated indices to avoid ambiguity, 
one obtains that 

mC M (x) = iK^t 0v x v - ^(iK) 2 J2 [Lo»,Uu // + 0(kV). (19) 

The covariant derivative can then be rewritten 

D„Q(x) = e**^ + r M (x))Q(x), (20) 

with the "connection" r M (x) now defined as 

T M (x) = k9,C M (x) = ~\^f [W, Uu] x v + (21) 

V 

It is easy to check that the first term of the expression is equal to the covariant 
curl, — (i/2)Fn U x u , which is non-zero due to the non-Abelian property of the gauge group 
50(4); it is also equal to the curvature tensor and in the present case is proportional to 
the gaussian curvature k 2 . 

When Q(x) = Q , the gradient term can be expressed at leading order in kI (with 
the Einstein summation convention again) as 

F cost (L) = \z f d 3 x | r M (x)Q | 2 

= iz K 4 (VQ )*(L^Qo) / d 3 xx v x v \ (22) 

o J(L) 

where a star denotes complex conjugation and where we have used the commutation 
relations of the 50(4) generators, [fjo^, Lo„ = ^L^. It is trivial to derive from 
the above equation that the free-energy cost grows as k 4 L 5 with some irrelevant 
pref actor [79^. 
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If, as argued above, frustration indeed generates a super-extensive free-energy cost 
opposing the extension of the ideal order, the associated effective interaction should 
be long-ranged: the integral of a coulombic interaction oc l/r D ~ 2 over a region of size 
L in dimension d precisely leads to an energy going as L d+2 . Such a Coulomb-like 
dependence is a very general feature of interactions between defects, be they vertices in 
superconductors, dislocations in crystals, or disclinations in bond-orientationally ordered 
phases jJ3|. Therefore, whereas the extension of the liquid LPS is more conveniently 
formulated in terms of a local order parameter, frustration is more conveniently 
formulated in terms of defect densities, akin to local "disorder" parameters. A full- 
fledged gauge field theory could possibly handle this duality (see above and Ref. |43j). 
However, a much simpler picture is provided by working with a single variable S(x) for 
all effective interactions. This leads to (possibly screened) Coulomb-frustrated models. 
In a continuum description, the associated free-energy functional is given by 

F[S] = J d*x{\z | ^S(x) | 2 +\r | S(x) | 2 +0(S 3 )} 

- l -K J d 3 x J d 3 yS(x)G(|x-y|)S(y), (23) 

where K > is a measure of the frustration strength. G(|x|) is the Coulomb 
interaction that behaves as l/|x| in 3d, possibly screened at distances much larger than 
an intrinsic frustration length £. (Typically, such a screening leads to a Yukawa pair 
potential, exp(— (|x|/£))/|x|.) In the above equation, we have been cavalier with the 
tensorial character of the variables and the interactions. 

One may as well consider a lattice version 

H=-jJ2 Si.S,- + y J2 S < G (* - x i) S i ( 24 ) 

(ij) hi 

which is more convenient for computer simulations. The models are chosen such 
that in the absence of frustration (K = 0), there is a continuous or weakly first-order 
transition to an ideally ordered phase at a temperature of order J. 

In the corrugated-space picture (33 | l?7 | ISO ] , one can interpret the above Coulomb- 
frustrated models as describing the competition between the extension of the ideal order 
in regions which look like the reference space, i. e., regions with an effective curvature 
that allows tiling by the LPS, and the repulsion between regions of alike curvature, which 
enforces the constraint that the average curvature must be zero (see also Ref.|%T]). 

Models similar to those described here have been used in quite different contexts. 
In diblock copolymers formed by two mutually incompatible polymer chains attached to 
each other, the repulsive short-range forces between the two types of components tends 
to induce phase separation of the melt, but total segregation is forbidden by the covalent 
bonds that link the subchains together [H21 IHSl IH1| • A microphase separation transition 
occurs instead at low enough temperature and the system then forms periodically 
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modulated phases, such as lamellar, hexagonal, and cubic phases. Similarly, self- 
assembly in water-oil-surfactant mixtures results from the competition between the 
tendency of oil and water to phase separate and the stoichiometric constraints generated 
by the presence of surfactant molecules, constraints that act as the electroneutrality 
condition in a system of charged particles [HHl IHEl EH IHHj- The same kind of physics 
also arises in a very different situation: stripe formation in doped antiferromagnets like 
cuprates has been ascribed to a frustrated electronic phase separation, by which a strong 
local tendency of the holes to phase separate into a hole-rich "metallic" phase and a 
hole-poor antiferromagnetic phase is prohibited by the long-range coulombic repulsion 
between the holes [89 . 

The behavior of Coulomb-frustrated models has been studied both through scaling 
arguments and by computer simulations. The results will be described in sections El and 

m 

5. Avoided criticality 

It seems a little distressing that no unique, unchallenged minimal theoretical model for 
frustration in liquids has been derived so far. However, the situation is not as bad as 
it looks, since there appears to be universal or to the least robust properties shared by 
the models described in the above section. The most important such property is that 
of "avoided criticality" [2011001 IHJj- This notion expresses the fact that the critical point 
which exists in the absence of frustration and separates a disordered (liquid) from an 
ideally ordered phase disappears as soon as an infinitesimal amount of frustration is 
introduced. 

Let us try to make this more explicit. Long-range ideal order is forbidden in 
the presence of frustration, but other types of order, which we generically describe 
as "defect ordered phases", are still possible. By defect-ordered phases, we mean phases 
in which the defects that break the ideal order themselves arrange in a periodic array. 
This includes: Frank-Kasper-like phases in materials with local icosahedral order, in 
which the disclination lines form a regular network; the Abrikosov flux lattice in type- 
II superconductors, in which the vortices form a lattice; and in Coulomb-frustrated 
systems, phases with modulated order, in which domain walls arrange in regular patterns 
leading for instance to lamellar or stripe phases. The phenomenon of avoided phase 
transition then means that the limit to zero frustration of the possible transitions to 
those defect ordered phases is at a temperature Tr>o{f rustration — > + ) significantly 
lower than Tc°\ that of the critical point in zero frustration. This discontinuity, which is 
illustrated in figure is a genuine non-perturbative effect due to the strong fluctuations 
generated by frustration. A similar phenomenon occurs if the transition at 

T (o) 

is 

weakly first-order [59] . (For a strong first-order transition, the transition may also be 
avoided, but the associated phenomenology is different [02] •) It is worth stressing that 
this property of avoided criticality has been overlooked in earlier studies of geometric 
frustration (201 EH]- 
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Frustration 



Figure 7. Schematic temperature-frustration diagram illustrating the property of 
avoided criticality (a genuine discontinuity when frustration goes to zero). 



5.1. Evidence for avoided criticality in the frustrated models 

We first discuss the cases where frustration is described either through a coupling to 
a non-Abelian gauge background or by an effective Coulomb interaction. Proving the 
existence of an avoided critical point then implies several steps. 

The first step is to show that the minimum of the interaction kernel A(|x — y|), or 
in Fourier space A(g), which appears in the quadratic part of the free energy functional, 
i. e., in 

F[S] = \J ^S(-q)A(g)S(q), (25) 

or a similar expression for the models defined in (Eq.(JHJ)), occurs at a non-zero wave 
vector with modulus go- m the Coulomb frustrated models (Eq. (j23j0 . the interaction 
kernel can be expressed as 

A(g) = q 2 + (AnK) /q 2 + constant, (26) 

where we have used that the Fourier transform of the Coulomb interaction is equal to 
47r/g 2 ; the minimum then occurs for go = {4txK) 1 / a (here and in the following we have 
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set for convenience Z = 1). In the screened Coulomb case, the above result is modified 
to 



with the minimum now at go = (4:TtK — f -2 ) 1 / 2 . For physical reasons (see section 
I4.3|) . one expects the coupling parameter K to be related to the inverse of the frustration 
length, K ~ so that the minimum is again at a nonzero wave- vector whenever 



For the models with a coupling to a gauge background, the demonstration is more 
involved. The interaction kernel in Fourier space is not only a function of q but also 
an operator or a matrix that should be diagonalized. It can be checked in the SO (4) 
theory of Nelson and coworkers (see Eqs. (fTJ) and Q) that all the eigenvalues \ l mm i(q) 
of the kernel, associated with a given representation of dimension (£ + l) 2 of 5*0(4), 
have a single minimum at a non-zero value qo, e mm ' that goes as k for small k (where n 
is the coupling strength associated with geometric frustration): see for instance figure 6 
of Ref . j2H] • In figure |H1 we also plot as another example the eigenvalues for the simpler 
case of 0(4) spins with 50(4) couplings in 3 dimensions |HH] . 

The second step consists in studying the models in the limit where the number 
of components N of the local order parameter is large. Exact analytic results can be 
obtained when N —>■ ooJSU] (this corresponds to the spherical limit for hard spins and 
to the gaussian variational or Hartree approximation for soft spins). A phase transition 
(of second-order in this case) occurs whenever 



where r is defined in Eqs. ()3l23|) and X(q) is the smallest of the eigenvalues of 
the interaction kernel. A similar expression can be derived for a lattice, hard-spin 
version [90 . It is easy to check that the integral on the right-hand side is dominated 
by wavevectors whose modulus q is near q Q ; in this domain, the denominator can be 
approximated by a quadratic term proportional to (q — q ) 2 whereas the integration 
measure is proportional to dqq^~ l \ so that when q ^ the integral strongly diverges 
in all dimensions d. Therefore, no ordering transition is possible at finite r, i. e., at 
non-zero temperature. 

When the model is considered on a lattice and not in continuum space, symmetry- 
breaking effects are generated that may induce a transition at a non-zero temperature. 
However, in dimension d < 4 this temperature is of the order of the magnitude of the 
symmetry-breaking terms, quite different and lower than the critical temperature T c 
in the unfrustrated system 90J. This leads to the phase diagram illustrated in figure [7| 
1/N corrections have also been computed and shown to leave unchanged the avoidance 
of the critical point [9"T]. 

As a third step, one can go beyond the large N limit (note however that the 
dimension of the main component of the order parameter for local icosahedral order in 



A(q) = q 2 + t~ 2 + (4nK)/(q 2 + t~ 2 ) + constant, 



(27) 



< w < 4. 




(28) 
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Figure 8. Eigenvalues of the quadratic part of the Hamiltonian describing 0(4) spins 
with non-Abelian 5*0(4) couplings in 3 dimensions in the continuum limit (see text). 
There are four eigenvalues, Ai = A 2 = fi — 6 + k 2 + q 2 , X^^ — fi — 6 + 2k 2 + q 2 T 
k\J k 2 + 4q 2 , where k is the frustration strength and [i is chosen so that the minimum 
value of the lowest eigenvalue is equal to zero. From top to bottom, A3,Ai = A2 and 
A4 divided by k 2 , plotted as a function of the dimensionless wave vector (q/n). The 
minimum in A3 occurs on a shell qo/n = V3/2. 

Nelson's 5*0(4) theory is equal to 169, quite a large value) and evaluate the magnitude 
of the fluctuations about the putative ordered ground state. This can be done in the 
usual and somewhat heuristic way by studying the quadratic fluctuations only. If the 
number of components N of the local order parameter is large enough, one always find 
"transverse" fluctuations whose size is given by the same integral as that on the right 
hand side of Eq. (|28| . Because of the strong divergence of the integral, the ground state 
is always unstable with respect to such transverse fluctuations: as a result, no ordering 
takes place at finite temperature as soon as frustration is present [231111 • In the Coulomb 
frustrated models, this is what occurs for N > 2. For N = 2, the divergence of the 
integral is only finite-temperature logarithmic, which suggests that a phase with quasi- 
long range order is possible |9 1.. . 

The case N = 1 (Ising), as more generally the case of all models with a discrete 
instead of a continuous symmetry, is also special. As first shown by Brazovskii[93 
and further supported by computer simulations |94| 193]. there is a first-order transition 
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to modulated phases; this first-order transition is induced by the strong fluctuations 
associated with frustration (indeed the mean-field approximation predicts a second-order 
phase transition) and its temperature Too goes continuously to Ti " 1 as frustration goes 
to zero. However, the limit of stability (spinodal) of the disordered phase is depressed 
to zero temperature and the critical point at remains the only finite-temperature 
critical point in the temperature-frustration diagram. The transition at T^o being first- 
order, it should be possible in principle to supercool the disordered (liquid) phase to 
lower temperature. This will be discussed in more detail in section UJ 

The above reasoning can be made more rigorous with the help of a generalized 
Mermin- Wagner inequality [SHI EE]; it shows that whenever the integral 



diverges, ordering at wave vector q is inhibited at finite temperature. 

Finally, some results are also available concerning the uniformly frustrated spin 
models. In the 2-dimensional version (Eqs. (JBJ), ©), there appears to be an avoided 
critical point, just like in the models described above. In this case, the defect-ordered 
phases are the Abrikosov vortex lattice and the associated hexatic vortex phase, and 
vanishing frustration means vanishing applied magnetic field. As shown by D. Fisher [75] 
on the basis of the continuum Coulomb gas description (Eq. ljllj) ) within the KTNHY 
theory of 2-dimensional melting (for a review, see [HZ]), the transitions to the defect- 
ordered phases occur, in the limit of vanishing frustration, at temperatures more than 
an order of magnitude lower than the Kosterlitz-Thouless transition [9$J EH] that takes 
place at zero frustration/zero magnetic field; this has also been confirmed by computer 
simulations [7T ] 172 ] HOOj . The resulting phase diagram, as illustrated by figure 1 in Ref. 
[75] . is then similar to that shown here in figure UJ 

Again, the physics of the transitions in the presence and in the absence of frustration 
is very different. The melting transitions describe the unbinding of the dislocations and 
disclinations that are present in the low temperature lattice or hexatic phases formed by 
the vortices induced by the applied magnetic field. On the other hand, the Kosterlitz- 
Thouless transition at T'J '' describes the unbinding of a neutral assembly of vortices 
of both positive and negative signs. The 3-dimensional version describing an extreme 
type-II superconductor in an applied field does not seem to have an avoided critical 
point jlUlj HU2j . but the case of an isotropic frustration discussed above in section IQ1 
has not been studied. 

5.2. Consequences of an avoided critical point 

In dimension d > 1 an Ising ferromagnet, when placed in a uniform external field, has 
of course some sort of avoided criticality: an infinitesimal field destroys the critical 
fluctuations beyond a given length scale. However, the phenomenon is trivial. The 
external magnetic field induces a non-zero magnetization at all temperatures, so that 
the system is always ferromagnetic and does not undergo any phase transition. 




(29) 
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Quite different is the situation in a ferromagnet with iV-component vector spins 
(and O(N) symmetry) in the presence of a quenched random field (with zero mean). For 
dimensions between 2 and 4, the pure system (no random field) possesses a critical point 
associated with the usual transition from paramagnetism to ferromagnet ism, whereas 
the random-field model, provided iV ^ 3, has no phase transition at all and stays 
paramagnetic at all temperatures |103*j : indeed, the lower critical dimension is d = 2 in 
the former case and d — 4 in the latter. Interesting, non-trivial scaling properties emerge 
in such a case for temperatures equal and less than Tc°\ the critical point of the pure 
system[103j. This is to such phenomenon that we refer when exploiting the property of 
avoided criticality shown by frustration-based models of glassforming liquids. 

On general grounds, one expects that the presence of an "avoided critical point" 
leads to several diverging length scales, according to whether one approaches the critical 
point at zero frustration by varying the temperature, or one approaches the critical point 
or the ideally ordered phase by decreasing the frustration to zero. A phenomeno logical 
scaling approach will be developed and discussed in the next section. 

6. Phenomenological scaling approach of glassforming liquids: 
frustration-limited domain theory 

The above developments suggest two main conclusions concerning the description of 
glassforming liquids 20J: 

1) Frustration naturally leads to collective (or cooperative) behavior on a 
mesoscopic scale, a feature that we stressed as being characteristic of the phenomenology 
(see sectional). The collective property comes from the phase growth (in which the liquid 
LPS spreads in space) induced by the proximity of an avoided critical point, whereas 
the limitation on the scale over such a growth can take place results from frustration 
that aborts the phase transition and leads to "avoided criticality". 

2) The relevant temperature about which one can organize a scaling description of 
the viscous slowing down and other collective properties of supercooled liquids is that of 
the avoided critical point in the associated unfrustrated system. This, of course, is only 
meaningful if the critical point is narrowly avoided, which implies that the frustration 
characterizing the liquid under consideration is small enough. Note also that in a liquid 
the temperature T^ marks a crossover, not a true transition, so that it is not expected 
to be sharply defined. 

6.1. Heuristic scaling arguments 

As a first step towards a scaling analysis, one can simply incorporate the physics of 
frustration-induced avoided criticality in the consideration of aborted nucleation of the 
ideal ordered phase in the liquid phase. At temperatures sufficiently below Tc°\ the free 
energy of a single domain of the ideal phase in a disordered liquid environment, when 
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its characteristic linear size is L, can be written as [201 ITU4*j 
F(L, T) = a(T)L e - 0(T)L 3 + s(T)L 5 . 



(30) 



The first two terms are those commonly encountered in classical nucleation theory: 
the first one represents the free-energy cost associated with building a domain of one 
phase in another; for simplicity, and in the absence of more specific information, we set 
6 = 2, which corresponds to a surface tension term associated with the interface between 
the two phases [TUB] . The second term is the bulk free-energy difference between the 
two phases and corresponds to a gain in free energy. If £o is the correlation length 
of the unfrustrated sytem, the surface tension a scales as ^ 2 and the ordering bulk 
free-energy density as £^~ 3 . The last term represents the strain free-energy resulting 
from the effect of frustration; as discussed in section 14. 3[ it is generically expected to 
grow super-extensively with system size, as L 5 in d = 3. The above equation assumes 
that the size L is large compared to the typical molecular length, so that a continuum, 
thermodynamic-like description can be used, but is small compared to the intrinsic 
characteristic scale of frustration, so that one does not have account to for the screening 
effect that modifies the super-extensive growth of the strain free energy at long distances. 

Because frustration in a liquid should be the same in any point of space, one 
actually expects, instead of the formation of a single domain, the breaking up of the 
liquid into a collection of domains, or, to use more pictorial terms, a " mosaic" 
or a "pathwork" |10fij of domains. Those domains are separated from each other by 
"interfaces", i. e., regions in which the ideal order is broken and a higher concentration 
of defects, such as disclination lines, is present. The size of the domains and their 
further growth when lowering the temperature is limited by frustration. Assuming 
that the inter-domain interactions are sufficiently weak to be treated by a mean-field 
approach, the free-energy density, $(L,T), for a domain of linear size L in a mosaic 
of frustration-limited domains (with randomly oriented local order parameters) can be 
directly derived from the free-energy of a single domain given above, namely f4*Hj. 



where a, 0, and s are only trivially renormalized with respect to their values for an 
isolated domain. 

The typical domain size L* is obtained by minimizing the free-energy density given 
above, and this leads to L*(T) ~ (a/s) 1 / 3 ; the strain coefficient s is, however, a priori 
unknown. The scaling analysis thus relies on (at least) two supermolecular structural 
lengths, the correlation length of the unfrustrated system, £o(^)> an d the typical domain 
size, L*(T). As one decreases the temperature below Tc°\ £o(^) decreases whereas 
L*(T) is expected to increase. Scaling only makes sense below a temperature 7\ < Tj ^ 
at which £o(^) an d L*(T) are comparable to each other and both larger than the typical 
molecular length. Such conditions may be only marginally satisfied in real liquids. 

It is important to stress that the mosaic of frustration-limited domains is an 
equilibrium feature (in the supercooled liquid phase it is of course metastable with 




(31) 
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respect to the actual crystal): the domains, as well as the overall pattern they form, 
are not frozen but continuously changing due to molecular motion; their statistical 
properties, however, depend only on the thermodynamic state of the liquid. 

In such a mosaic of domains, a relaxation and flow must involve restructuring of 
domains. In the simplest scaling picture, the typical free-energy barrier to be overcome 
for such a restructuring, AE*(T), scales as a(T)L*(T) 2 , i. e., involves motion, creation, 
or destruction of a domain wall. A more detailed argument for such a description has 
been given by Stillinger in terms of a "tear and repair" mechanism[106j. Scaling about 
a narrowly avoided critical point then implies that both cr(T) oc £o(T)~ 2 and L*(T) have 
a power- law behavior in (T^ — T) for temperatures less than Tc°\ (Recall that at and 
above there are no frustration-limited domains.) As a result, the typical free-energy 
barrier for a relaxation scales as 



In a naive, mean-field-like picture, o~(T) and L*(T) vanish as one approaches T c 
from below in an analytic way, o~{T) oc (T^ — T) and L*(T) oc (T^ — T), so that one 
obtains ip = 3. In a more elaborate, but still heuristic [l()5j description , it has been found 
that L*{T) goes as K~ 1 / 2 ^q{T)~ 1 , where K is the relative amplitude of the frustrating 
interaction (20] (see Eqs. (j23| l2"l*|) ): the typical domain size then grows as 
where v is the correlation length exponent of the unfrustrated system, as temperature 
decreases; it also increases as frustration decreases and diverges in the limit of zero 
frustration as expected in a scenario of avoided criticality. Taking v m 2/3, which is 
characteristic of ordinary critical phenomena in the absence of quenched disorder in 3 
dimensions, one arrives in this case at ip ~ 8/3. 

The above estimates concern the collective contribution to the free-energy barrier 
for a relaxation, that due to the frustration-limited domains which vanish around and 
above Tc°\ One must also account for the molecular contribution which characterizes 
"ordinary" liquid behavior (typically, above the melting temperature T m ). The simplest 
description that provides a reasonable fit to the experimental data in the ordinary liquid 
phase without introducing any spurious singularities is a plain Arrhenius formula, 



where and are species-specific molecular constants characteristic of the high- 
temperature liquid. 

The main prediction of this scaling approach is thus that the a-relaxation time and 
the viscosity have an activated-like dependence on temperature, with a crossover from 
molecular, ordinary liquid behavior to collective, domain-dominated behavior around 
T* » Tc°\ a crossover that leads to a change from Arrhenius- like to superArrhenius 
T-dependence, 



AE*(T) ~ (Tp -Tf. 



(32) 




(33) 




(34) 
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with E(T) = + AE*(T). The superArrhenius contribution to the activation 
free energy is described by a universal power law, 

AE*(T) = forT>T c {0) , 

= Bk B T^(l-^f for T<Ti°\ (35) 

This prediction fits well the extant experimental data on glassforming liquids with 

t/> in the range from 7/3 to 3 (the best overall fit being for tp = 8/3)1201 El) a range 

compatible with the above estimates. This is shown in figure El In the fitting procedure, 
T (o) 

is an adjustable parameter: for virtually all liquids studied, its value is found slightly 
above the melting point T m ; this indeed makes sense since one expects the stability of 
the ideal ordered phase in the unfrustrated system at to be higher than that of the 
actual crystal in the frustrated liquid jS3 EH] ■ 

The parameter B in Eq. f)35j) is a measure of the departure from Arrhenius behavior, 
hence of "fragility" . From the scaling analysis, we obtain that B increases as frustration 
decreases (it goes as K^ 1 for Coulomb frustrated models). One can thus associate a large 
fragility with a small frustration. This is compatible with the view developed in sectional 
that atomic glassformers made of spherically symmetric particles are not fragile, or in 
the standard terminology are rather "strong", despite the weakness of the inter-particle 
interactions: frustration in those systems characterized by local icosahedral order is 
indeed strong. 

6.2. Further dynamic and static scaling analysis 

It is tempting to push further the scaling analysis in order to address other aspects of 
the phenomenology of supercooled liquids. 

The equilibrium distribution of domain sizes in the mosaic, p(L, T), at temperatures 
below Ti < Ti 0) can be obtained from p{L,T) oc exp [-($(£, T) - <$>*(T))L 3 /k B T)\, 
where <3>(L,T) is the free-energy density given in Eq. (j31j) and $*(T) its value at the 
minimum. It can be cast in a scaled form as |481 1107] 

p{L,T) ex exp{- 7 (T)[K(A) 2 - 3 -(±f + \(^)% (36) 

where the entire temperature dependence is contained in the parameter j(T) oc 
k is a number of order 1 introduced to allow for small domain- 
shape effects at a mean-field level 1 1 OTI (it should not be confused with the frustration 
parameter k used in sections El and HJ. 

By generalizing the simple dynamic scaling arguments presented in the above 
subsection, one can also evaluate the collective contribution to the activation free energy, 
AE(L,T), for a domain of size L; this latter takes the following scaled form [4"Hl I107j : 

AE } L ' T) = & 7 (T)[(^) 2 - m(^)% (37) 
k B T n ny L* J K L* J J v ; 

where b and m are numbers of order 1, with bm < 1/2. The last term in L 5 
has been included to account for the fact that building or moving a domain wall costs 



The frustration-based approach of supercooled liquids and the glass transition 



31 



□ 14 - 



?2- 



0. lO- 



Cl OS - 



X 

1 

t- C C4 - 

02- 



D UD 




Figure 9. Scaling of the collective component of the effective activation (free) 
energy, AE* (T) , versus reduced temperature (T* — T) /T* measured from the crossover 
temperature T* (or equivalently, the avoided critical temperature ) . The solid line 
is the 8/3 power law discussed in the text. Reprint from Ref.0]. 



a free energy of order aL 2 but leads to a reduction in strain free energy. The full 
effective activation free energy is then obtained by adding the molecular, Arrhenius-like 
contribution, E^/kBT, to the above expression. 

Below the temperature T\ < at which the liquid breaks up into supermolecular 
domains, the dynamics are naturally described as heterogeneous. For instance, the 
nonexponential character of the dielectric relaxation can be simply modeled by assuming 
that molecular reorientations are completed within a single domain, in which they are 
dynamically coupled to the order variable, and by assuming that the relaxation within 
each domain is exponential in timej^] . (Of course, there may be systems such as polymers 
for which the bare "molecular" relaxation is already non-exponential; the collective 
behavior associated with the mosaic of domains then increases the non-exponential 
character.) This leads to the following expression for the normalized dielectric relaxation 
function UHllMj: 



f a (t) = / dLL p(L, T) exp exp — ^ '- 

Jo I Too V Kb J- J 



(3f 



with p(L,T) and AE(L,T) given in Eqs.(|3T)]) and (J37J). The above expression 
compare well with experimental data on glassforming liquids[107j. This is illustrated in 
figure [TUl where the imaginary part of the frequency-dependent dielectric susceptibility 
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Figure 10. Frequency-dependent dielectric susceptibility of the fragile glassformer 
m-fluoroanilin for several temperatures on a log- log plot: comparison between 
experimental data (T g = \72K) and frustration- limited domain scaling theory (data 
and fits by C. Alba-Simionesco, see also Ref. |108j '). The parameters used for the fit 
(see text) are T* = T c (0) = 258K, E x = 1587K, = 10" 12 s, « = 0.62, m = 0.29, 
b = 1.04, B = 364. The exposant ip is fixed to 8/3. 



of the fragile glassformer a m-fluoroanilin is compared to the Fourier transform 
of df a (T)/dt obtained from Eqs.(|5fi|). (j3*Tj) and (J55j) : the main features of the a 
relaxation are well reproduced over a large frequency range (up to 13 decades) and 
a large temperature domain [TUH] . (The approximate scaling plot proposed by Nagel 
and coworkers 109 is also well reproduced |48| ll()7j .) One of course recovers for the 
temperature dependence of the peak frequency the power-law behavior described in 
Eq. (j33j) . The number of adjustable parameters strongly restricts the predictive power 
of the scaling approach, but it should be stressed that those parameters are independent 
of temperature. From those parameters one can estimate the typical domain size, which 
is found to be around 5 to 10 molecular diameters at T g , an estimate that is compatible 
with experimentally determined sizes (see section EJ); with additional input [TTU] . one 
can also describe the decoupling between translational diffusion and reorientational 
relaxation |107| . 

A few additional points are worth stressing: (i) in this approach, the heterogeneous 
character of the a relaxation in the supercooled liquid regime is a consequence of a 
structural property, namely, the break-up of the liquid into a mosaic of frustration- 
limited domains; this is to be contrasted with other descriptions that focus on purely 
kinetic mechanisms (see the review article [TTT] ). (ii) The superArrhenius growth of 
the relaxation time with decreasing temperature is accompanied by only only a modest 
stretching of the relaxation functions (compare for instance with the quite different 
behavior of magnetic systems in the presence of a random field |112 l ; in the present 
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scaling description, this results from the decreasing polydispersity of the domains as 
temperature decreases. 

Finally, scaling about an avoided critical point gives a qualitative account of the 
rapid decrease of the excess (or " configurational" ) entropy observed in supercooled 
liquids as one approaches the glass transition. Within this approach, there are indeed 
two sources of entropy decrease below Tc^ as one lowers the temperature: the first one is 
the increase of ordering that should take place inside all the domains |113j . the second is 
associated with the growth of the typical domain size and the resulting decreasing weight 
of the interfaces with high concentration of defects. This latter effect will be further 
discussed in section |H1 in connection with the description of glassforming liquids in terms 
of the complexity of their free-energy landscape and the large number of metastable 
states. Just as the a-relaxation time and the viscosity rapidly increase with decreasing 
temperature with no need for a low-temperature singularity, the frustration-limited 
domain approach predicts a rapid decrease of the entropy with no requirement for a 
low-temperature Kauzman catastrophe. 

The above phenomenological scaling analysis illustrates how the concept of 
frustration coupled with the property of avoided criticality can provide a framework 
for describing the phenomenology of glassforming liquids. However, whether or not 
the assumptions that underlie the scaling analysis evolve rigorously from the statistical- 
mechanical models discussed in section EJremains a fundamental theoretical uncertainty. 
This point has not yet been resolved, but some progress has been achieved by means of 
computer-simulation studies of simple frustrated models. This is what we discuss next. 

7. Computer simulations of simple systems: Coulomb frustrated lattice 



The Coulomb frustrated spin models described in section 14.31 seem to be the simplest 
systems compatible with the frustration-based approach of glass formation. Setting 
aside the issue of their microscopic derivation, it is possible to check the relevance of the 
scenario presented above by investigating numerically such models. Extensive Monte- 
Carlo simulation studies [551 II 14[ ITTo] have been performed for models on a 3-dimensional 
cubic lattice with Hamiltonian (see also Eq.(J2U) ) 



for a variety of spin variables, mostly with discrete symmetry: Ising model, S; = 1, 
g-state clock model, Sj = (cos(27r#;/g), sin(27r#j/g)), with 6{ the orientation of the planar 
spin and q = 5 and q — 11. We also include here new results obtained for the XY model, 
i. e., for continuous planar spins. 

In the absence of frustration (Q = 0), the models have a critical point at a 
temperature Tc°\ below which they are ferromagnetic. The long-range, coulombic 
interaction requires that the total magnetization of the system be zero in order to 



models 




(39) 



The frustration-based approach of supercooled liquids and the glass transition 



34 





Figure 11. Spin-spin correlation function C(t) versus ln(t) for several temperatures 
both above and below (in each panel, curves from left to right are for decreasing 
temperatures): 5-state clock model for Q = 0.1 (a) and Q = 0.00625 (b); 11-state 
clock model for Q = 0.001 (c); and XY model for Q = 0.005 (d). 
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ensure a proper thermodynamic limit. Therefore, long-range ferromagnetic order is 
prohibited at all T's for any nonzero value of the frustration parameter Q/J. (In 
what follows, we take J = 1.) The models studied (XY and discrete orientations) 
display only a "weak" form of avoided criticality. They can still form ordered or quasi- 
ordered phases characterized by modulated patterns (e.g., lamellar phases) below a 
temperature T DO (Q) that goes non-analytically but continuously to as Q goes to 
zero. However, the transition between paramagnetic and modulated phases is first- 
order [HUES]- Following Brazovskii[93 , this result can be interpreted on the basis of the 
self-consistent Hartree approximation which predicts the occurrence of a fluctuation- 
induced first-order transition, a transition with no nearby low-T spinodal. As a result, 
T (o) 

remains the only nearby (and avoided) critical point in the frustration-temperature 
diagram. All the simulations have been made in the disordered (paramagnetic) phase. 
Since the transition to the modulated phases at Tdo(Q) is first-order, one could in 
principle supercool the paramagnetic phase to lower temperatures; however, it has 
been found that the lattice sizes accessible in practice are too small to allow for a 
proper supercooling below Tnn \114\ I115j . This, unfortunately, restricts the domain of 
temperature over which one can study the slowing down of the dynamics. (Out-of- 
equilibrium dynamics and aging phenomena have not been considered.) 

The relaxation to equilibrium of the frustrated systems has been studied via the 
Monte Carlo algorithm (Metropolis rule), time being the number of sweeps per spin. The 
dynamical quantity that has been monitored is the equilibrium spin-spin self correlation 
function, C(t) = (1/N) J2i < Si(£')Si(£' + 1) >, where the bracket denotes the thermal 
average and N is the total number of lattice sites; cubic lattices of size 16 3 to 22 3 with 
periodic boundary conditions have been used and the Coulomb interaction has been 
handled via Ewald sums |1 14*1 1115j . 

The evolution of C(t) with T is illustrated in figure fTTI for the 5-state clock (Q — 0.1 
and Q = 0.00625), the Ising (Q = 0.001), and the XY (Q = 0.005) models. At the 
highest frustration, illustrated by figure ITTk. the decay of C(t) appears to proceed in a 
single step at all temperatures. For lower frustrations, a 2— step decay develops as T is 
lowered (see figure [TTJo-d). At high temperature, typically above the critical point of the 
corresponding unfrustrated model (T^ ~ 2.1, 4.51, 2.32 for the 5— state clock, the Ising, 
and the XY models, respectively), and for all frustrations, the whole time dependence of 
C(t) is well fitted by a simple exponential, exp(— t/r (T)). At low temperature, below 
Tc°\ it is impossible to describe the entire decay of C(t) by a single exponential; the 
emerging second step of the relaxation can be described by a stretched exponential, 
exp(— (t / TKwwiT)) 13 ^) , while the first step can still be fitted by a simple exponential, 
expf— t/Tn(T)) frn\ I115j . This non-exponential behavior and emergence of a two-step 
decay, both of which become more pronounced as T is decreased, are typical of fragile 
supercooled liquids. One may however notice that the timescale separation between the 
two relaxation steps that can be achieved in the simulations is not sufficient to observe 
the development of a true plateau at intermediate times. 
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Figure 12. Arrhenius plot of r a (T) for the 5-state clock model and various 
frustrations (from left to right in (a) and from right to left in (b), Q = 
0.00625,0.00125,0.025,0.05,0.1). In (b) the curves are plotted versus the reduced 
inverse temperature T s (ln(r) = 5.9) /T, where T g (ln(r) = 5.9) is the temperature at 
which the maximum relaxation time, ln(r a ) = 5.9, is attained. 
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The slowing down of the relaxation as one lowers the temperature is visible in 
figure The temperature dependence of the relaxation times characteristic of the 
short- and long-time behaviors, tq(T) and r a {T) respectively, has been investigated in 
Refs. |114] I115j . Over the range of temperature studied, the dependence of r (T) is 
essentially Arrhenius-like, t (T) ~ r 0iO o exp(£' 00 /T), whereas the dependence of r a (T) 
shows a marked deviation from Arrhenius behavior below some crossover temperature 
in the vicinity of Tc . 

The crossover from Arrhenius to super- Arrhenius behavior of r a is shown in figure 
IT21 for the 5— state clock model for several different frustrations. A similar trend is 
observed for the Ising and XY models. Several points are worth noting: (i) such a 
super- Arrhenius behavior is typical of the viscous slowing down of fragile glassforming 
liquids, the more fragile a liquid the more pronounced the super-Arrhenius character; 

(ii) the crossover occurs in the vicinity of the critical point of the unfrustrated system; 

(iii) the departure from Arrhenius behavior, i.e., the "fragility", becomes more marked 
as frustration decreases. 

As shown in Refs. |114j and |115j , the temperature dependence of the a-relaxation 
time is compatible with the behavior predicted by the frustration-limited domain 
approach (see Eq. (J35J)). However, because of the limited domain of temperature, 
hence the limited range of relaxation times, one should not put too much emphasis 
on the fit: for instance, the data are also compatible with a Vogel-Fulcher formula, 
r a = r o ex P(DT /(T — Tn)) |116j : on the other hand, they cannot be fitted by a power- 
law dependence. 

Although no direct determination of a domain or heterogeneity length has been 
attempted in those studies, interesting insight has been provided by investigating the 
finite-size effect on the relaxation time. As illustrated in figure IT3*k. the behavior is quite 
different than that of standard critical slowing down (shown here for the unfrustrated 
system, figure fTTTb) : in the presence of frustration, the slowing down of the relaxation 
does not come with the rapid growth of an associated length, which is compatible with 
the experimental observations on glassforming liquids. 

The computer simulation studies of Coulomb frustrated spin models thus provide 
some evidence that such models display many of the main characteristic features of the 
phenomenology of fragile glassforming liquids. It is however fair to say (i) that the 
analogy with supercooled liquids is limited by the range of relaxation times accessible in 
simulations (due to the presence in those models of a first-order transition to modulated, 
defect-ordered phases) and (ii) that the mechanisms by which frustration induces a 
slowing down of the relaxation remain at the level of speculations. 

8. Connection to other approaches of the glass transition: complex 
free-energy landscape and kinetic constraints 

So far, we have stressed how frustration provides a physical mechanism for generating 
cooperative behavior on a mesoscopic scale, a feature that seems to be at the very root 
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Figure 13. Arrhcnius plot of r a for the 5— state clock model and for lattice sizes 
L = 5,10,15,20 (from right to left): (a) Q — 0.00625 (no difference is seen between 
the 3 larger sizes) and (b) Q — (standard critical slowing down). 



of the phenomenology of fragile glassforming liquids. We have focused on a generic 
property induced by frustration, namely avoided criticality, and on the possibility of 
developing a scaling approach of supercooled liquids exploiting this property. However, 
starting from the same statistical mechanical description based on frustration, other 
routes are possible. It is fruitful to make the connection with other, maybe more 
common, lines of thought concerning the emergence of glassy dynamics: the "landscape 
paradigm" , which relates the slowing down of the dynamics to the presence of a large 
number of long-lived metastable states, and the approaches in terms of "kinetically 
constrained models" , in which glassiness results from constraints on the effective kinetics 
of the system. Whether or not one description is better than the others and to what 
degree the various approaches are compatible one with another are still, we believe, 
open questions. 

8.1. Long-lived metastable states and complex free- energy landscape 

On general grounds, one expects frustration to give rise to a large degeneracy of low- 
lying energy configurations \117\ 126 j . A first insight into the presence of long-lived 
metastable states in supercooled liquids can be gained from a heuristic argument that 
follows from the phenomenological scaling approach discussed in section |U1 If indeed 
below some temperature a glassforming liquid splits into a mosaic of frustration- 
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limited domains, if such domains or small groups of domains are characterized by a 
limited number of long-lived metastable states whose exploration via thermal activation 
gives rise to the a relaxation, and if those domains are only weakly interacting (see 
section I6.2J) . then the system as a whole is characterized by an exponentially large 
(in system size) number of long-lived metastable states: those states are obtained by 
combining the metastable states of the roughly independent mesoscopic regions, which 
does lead to an exponential total number. Since the typical domain size is predicted 
to increase when either T decreases (at constant frustration) or frustration decreases 
(at constant T), the number of metastable states is expected to decrease when T or 
frustration decreases. The logarithm of the number of long-lived metastable states can 
be used to define a " configurational entropy". A non-zero " configurational entropy" 
thus emerges below the crossover temperature associated with the avoided critical 
point and decreases as temperature further decreases. 

How can one check the above property? The difficulty is twofold. First, a 
precise definition of a "long-lived metastable state" must involve a time scale and, as 
a consequence, cannot a priori be treated within a purely static calculation; secondly, 
one has to find an operational way to compute the number of such metastable states. 
In the absence of a better solution jl!8j . a way to to handle the first aspect is to assume 
that the long-lived metastable states, if indeed present in the system, have an infinite 
lifetime: such an approximation, which neglects the thermally activated processes that 
allow the system to escape from the metastable states and are the source of the a 
relaxation, has an intrinsic mean-field nature. The other aspect of the problem has 
found an elegant solution proposed by Monasson |119j : it relies on introducing a weak 
coupling to a pinning field that prevents the system from sampling all the metastable 
states (hence breaking ergodicity) and handling the average over the pinning field by 
means of the replica method, much used in theories of spin glasses and other systems 
in the presence of quenched disorder fl Q I12()j . 

When applied to the models introduced in section |S] to describe frustration in 
liquids, this "replica mean-field approach" proceeds as follows [T2*T] I122j . Consider a 
frustrated model described (in the continuum limit) by the Hamiltonian or free-energy 
functional F[x] (x( x ) stands for Q(x) in Eq.fjHJ) , S(x) in Eq. (j23|) . etc) and add a 
weak coupling to a pinning or "ergodicity-breaking" field cr(x). The associated modified 
partition function is given by 



Z x [a] = J V X exp -^F[ X ] - ~ j "fx (<r(x) - X (x)) 5 



Vxexp 
exp ( -(5 f\ [cr 



(40) 



with (3 = l/{k B T) and A > 0. If F[x] possesses low-lying minima, Z\[o\ will 
be large at low temperature when cr(x) is close to a configuration of the field x( x ) 
corresponding to one of those minima. A weighted average of the free energy in the 



The frustration-based approach of supercooled liquids and the glass transition 40 



various metastable configurations is then obtained as 



T ^) v °™j^^v (41) 

When A is equal to zero, T\ [a] reduces to the thermodynamic free energy of 
the system, T = — fcsTln(J Vxexp [— (3F[x]\) and so does F\=o- The behavior as 
A —>■ + may however be subtle: if the number of metastable states (low-lying minima) 
is exponentially large in system size, the approach A — ► + is non-perturbative and a 
discontinuity F\^o+ —F\=o is present; the magnitude of the discontinuity is precisely, up 
to a factor T, the "configurational entropy", defined as the logarithm of the number of 
metastable states. Again, this can only happen in a mean-field description, the system 
in the limit A — > + being stuck forever in one of the exponentially many metastable 
states. 

The replica trick can now be used to compute the average over the pinning field in 
Eq-dHJ. Defining 

prtm) = -— lnZ(m) = -— lim x ^+ [ VoZ x [a] m , (42) 
m m J 



m=l 
m=l t 



one obtains after straightforward manipulations that JF A ^ + = d(mT(m))/dm 
and the configurational entropy, S c = (3 (j r \_> + — T \ is equal to f3dJ 7 (m)/dm 
whereas the replicated partition function Z{m) can be written as 

/m m , m „ 

l[Vx a expi-pYs F {Xa} + ^ E / ^Xa(x)x6(x)]. (43) 

a=l a=l a,b=l 

The actual computation of Z(m) requires of course the introduction of 
approximations, mean-field-like approximations as stressed above. The first such 
calculation on a frustrated model has been performed by Schmalian and Wolvnes |121j 
on the Coulomb frustrated model (Eq. (|23Jl ) with an approximation known as the 
self-consistent screening approximation fSCSA) [l23| I124j They find that below a 
temperature an extensive configurational entropy, i. e., an exponentially large 
number of metastable states appear; the entropy decreases as T decreases and becomes 
negative below a temperature Tk- 

These results call for several comments: (1) they confirm the heuristic picture 
described above and, as expected, T4 is always found below Tc°\ the critical point of the 
model without frustration; one also obtains that the configurational entropy decreases 
at constant temperature when frustration decreases. (2) The scenario is similar to that 
characterizing generalized spin glasses, such as the p-spin of the Potts glass models, 
in the mean- field limit |1 20] . In those models, the upper temperature is associated 
with a dynamic singularity described by the ideal mode-coupling approach whereas 
the lower temperature Tk corresponds to a bona fide thermodynamic transition to a 
spin-glass phase. It can actually be shown that the temperature Ta in the Coulomb 
frustrated model is also associated with a dynamic singularity similar to that of the 
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ideal mode-coupling description |125| . One must however recall that this occurs within 
mean-field approximations and that such a dynamic singularity should vanish in an 
exact description. The analogy with generalized mean-field spin glasses is the basis of 
the random first-order transition theory of glass formation developed by Wolynes and 
his coworkers [T7| ITH| 11261 \1'27\ I128j . We shall come back to this theory below. 

The same property, namely the appearance of a complex free-energy landscape 
with an exponentially large number of metastable states has been found in other 
systems related to the frustration-based description of glassforming liquids, systems 
in which the minimum of the interaction kernel in Fourier space occurs at a non-zero 
wave vector g n [130. 18T | 11311 1X32"] . As discussed in section these systems all show the 
property of avoided criticality. However, most of the studied models are akin to the Ising 
version (i. e., with a discrete Z2 symmetry) and display the "weak" form of avoided 
criticality: a fluctuation-induced first-order transition. We have already mentioned that 
this transition, whose location is non-universal and model-specific, is hard to bypass by 
supercooling, which limits the accessible domain of temperature in computer simulation 
studies. As a result, the conclusions concerning the glassy behavior of the models may 
not be as clear-cut as one would like and details concerning the parameter range studied 
may become important: this may explain the apparently conflicting results obtained 
for the 3-dimensional Brazovskii-like model [T3*3*l 1125} Il34[ I135j . It seems anyhow that 
models with a continuous symmetry, like the Coulomb frustrated O(N) model, also 
follow the same scenario with the emergence of a complex free-energy landscape below 
the (truly) avoided critical temperature T^ [136J. 

One may wonder at this point how reliable are calculations based on approximations 
which are to a large extent uncontrolled? There is of course no definite answer to this 
question, but one can at least try to check the robustness of the scenario: what happens 
when improved, or on the contrary worse, approximations are used? The approximations 
being mean-field in character, they may overestimate, as is often the case, the presence 
of transitions and singularities: are there models for which the approximations do not 
predict the appearance of a complex free-energy landscape? The available answers to 
those questions tend to support the validity of the scenario. An exponential number 
of metastable states is also predicted when using an approximation developed for 
studying strongly interacting systems with predominantly local correlations |137j . the 
" dynamic mean- field approximation" [138 ; on the other hand, a number of perturbative 
treatments cruder than the SCSA mentioned above do not find a complex free-energy 
landscape |125j . Perhaps more interestingly, when the replica mean- field approach and 
the SCSA are applied to the 3-dimensional model in the absence of frustration, it does 
not find the signature of an exponentially large number of metastable states [119J, as 
indeed expected for a pure 4 theory. 

To conclude this section, it is worth summarizing the main points of the random 
first-order transition theory of glassforming liquids [T71 HHJ 11261 11^71 11281 EH]- m this 
approach, the mean-field scenario outlined above is modified to include thermally 
activated escape from metastable states. The dynamic mode-coupling-like singularity is 
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then smeared out to become a crossover to an activation dominated regime associated 
with the evolution of the system on the complex free-energy landscape. In real space, 
the liquid is described as a mosaic of cells and activation events as entropic droplets, the 
driving force for nuleating such droplets being provided by the non-zero configurational 
entropy. Super- Arrhenius temperature dependence of the a-relaxation time in the form 
of a Vogel-Fulcher expression is predicted as one approaches the ideal glass transition 
Tx', at this point the configurational entropy vanishes and it can therefore be identified 
to the Kauzman temperature [TT]. This approach successfully describes many aspects of 
the phenomenology of supercooled liquids. 

8.2. Defects and effective kinetic constraints 

Topological defects play an important role in frustration-based approaches of 
glassforming liquids. In the description of frustrated icosahedral order in 3-dimensional 
space for instance, the most relevant defects are disclinations, i. e., rotational line defects 
(see section |3J). We have also illustrated in section 0] how the statistical mechanical 
treatment of frustrated systems may become more tractable by going from the original 
representation in terms of local order parameter to a dual description in terms of 
topological defects. 

It must be emphasized, however, that the simplicity one could gain via such a 
transformation for studying the structural and thermodynamic properties may have a 
penalty associated with an increased complexity in the effective equation of motion 
for the defects. In general, one would have to deal with the kinetics of defect lines 
in a partially ordered but non-crystalline medium. The slow dynamics of defect lines, 
such as vortex lines in 3-dimensional extreme type-II superconductors, is controlled by 
crossing, cutting and reconnection, all processes that are thermally activated and may 
lead to an effective freezing in an entangled "glassy" state on the experimental time 
scale |139j . Entanglement effects are expected to be especially severe for disclination 
lines in a frustrated icosahedral medium: in such a case, the underlying theory has 
a non-Abelian SO (4) symmetry and the resulting combination rules for the defect 
(disclination) lines has been predicted to lead to non-trivial kinetic constraints [3~H I26j. 
The phenomenological scaling approach reviewed in section |H1 side steps the issue by 
merely considering quasi-ideally ordered regions (domains) with a low density of defects 
separated by regions (interfaces) with a high density of defects, without specifying the 
nature of the defects nor their kinetics. 

To shed light on the subtle aspects of the duality between the description in terms 
of local order parameter and that in terms of defects, it is interesting to consider 
simple models which do not belong to the class of frustrated systems discussed in 
section 0D They consist in Ising spin models on 2-dimensional lattices, with short-ranged 
ferromagnetic multi-spin, "plaquette"-type, interactions. In the original formulation 
involving the spin variables, the system is strongly correlated |14()^ IT4~T] (with trivial pair 
correlations but non-trivial higher-order ones), and the dynamics is described via local 
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kinetic rules involving 1-spin flips. In the dual version in terms of topological defects, 
the thermodynamics reduces to that of an ideal gas, but the effective dynamics of the 
defects is now strongly constrained |14()[ I14l| fl42j). 

The presence of kinetic constraints, even in the absence of non-trivial static 
properties, has been shown to induce a slowing down of the relaxation and, as 
a result, formation of glassy states on some observation time scale. There is an 
abundant amount of work on such "kinetically constrained models" (for a review, see 
Ref. which have been, most notably in recent years [T5| I143[ ITlj . advocated 

as an alternative to thermodynamic approaches for explaining the phenomenology of 
supercooled liquids. In such descriptions, glassy behavior is attributed to the emergence 
of dynamic heterogeneities, with a typical length scale that grows in space-time, but 
not in space, when temperature is decreased or density increased, the models reproduce 
many features of glassforming liquids, with however some shortcomings. In all studies 
so far carried out, the defects are point-like objects; defects becoming more and more 
dilute as one approaches the glass transition region, their contribution to the entropy 
and heat capacity of the liquid (in excess to those of the crystal) is then completely 
negligible, and the empirical correlation between slowing down of the relaxation and 
decrease of the "conflgurational" (excess) entropy observed in liquids approaching their 
glass transition cannot be described |144[ I145j . In addition, the origin of the kinetic 
constraints is quite elusive. 

It is then tempting to speculate, after Palmer et al. |146j that "it is possible to 
see static frustration over many length scales as the underlying cause of dynamical 
constraints over many time scales" . Defects could be line-like objects instead of point- 
like ones, which could help keep the connection between dynamics and thermodynamics. 
However, no such description has been attempted so far. 

9. Conclusion 

In this article, we have reviewed the frustration-based approach of supercooled liquids 
and the glass transition. Frustration in liquids, namely a ubiquitous incompatibility 
between the spatial extension of the locally preferred liquid structure and the tiling of 
the whole space, provides a physical mechanism to generate what we have argued to be 
the distinctive characteristic of the viscous slowing down leading to glass formation: the 
emergence below some crossover temperature of a cooperative behavior whose extension, 
however, remains limited to a mesoscopic length scale. Cooperativity is associated with 
the (collective) extension of the locally preferred structure whereas the limitation is a 
result of frustration. Glassiness is then self-generated as one decreases temperature. 

Despite the absence of a well established minimal theoretical model of frustration 
in glassforming liquids, progress has been made over the last decades in developing 
a statistical mechanical framework. Models have been built, based on such ideas 
as a coupling to a non-Abelian background, a uniform frustration, or a competition 
between effective interactions. Interestingly, several generic properties characterize such 



The frustration-based approach of supercooled liquids and the glass transition 44 

frustrated models, largely independently of their detailed structure. Those properties 
are: (1) avoided criticality, a genuine non-perturbative phenomenon induced by 
frustration, (2) the emergence of a complex free-energy landscape with an exponentially 
large number of "long-lived" metastable states, and (3) the presence of topological 
defects. 

The existence of the above three features open the way to several theoretical 
routes to explain the phenomenology of glassforming liquids within a frustration- 
based description. The presence of a narrowly avoided critical point in the diagram 
frustration-temperature has been used to develop a phenomenological scaling analysis of 
supercooled liquids, the frustration-limited domain theory. The temperature associated 
with the avoided critical point marks a crossover between "ordinary" liquid behavior 
and "anomalous" (cooperative, heterogeneous, super- Arrhenius) liquid behavior, and 
fragility is predicted to be inversely related to frustration. Along a different line 
of thought, the emergence of a complex free-energy landscape is at the basis of the 
random first-order transition theory in which (super-Arrhenius) activated relaxation 
is associated with the configurational entropy resulting from the exponentially large 
number of metastable states. Finally, the existence of topological defects (primarily, line 
defects in 3-dimensional geometrically frustrated media), which can form an entangled 
state and whose dynamics may be strongly constrained, is suggestive of a possible 
treatment by means of generalized kinetically constrained models. Whether or not 
these theoretical frameworks are compatible and which of them is the most efficient to 
describe the physics of glass formation are open questions. 

The phenomenological scaling approach about the putative avoided critical point 
leads to a good qualitative and even, at the expense of course of introducing adjustable 
parameters, quantitative description of the salient phenomena in glassforming liquids. 
However, its theoretical foundations, besides the property of avoided criticality, remain 
shaky. Support has been provided by computer simulation studies of Coulomb frustrated 
models, but no definite conclusion should be made. Needless to say that in spite of its 
appealing features and of the results derived so far, much still remains to be done to 
make the frustration-based approach of supercooled liquids an operational, predictive, 
and microscopically founded theory of glass formation. 
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